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Abstract. I discuss some aspects of the Disc Instability Model applied to outbursts in vari- 
ous types of binary systems. I lament the general lack of interest in the subject. 
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1. Introduction 

Disc outbursts are observed in all kind of bi- 
nary systems. The best observed are outbursts 
of dwarf-nova stars where the disc fed by a red- 
dwarf companion surrounds a white dwarf pri- 
mary star Initially, there was some confusion 
about the dwarf-nova eruption site but Smak 
(1971) clearly identified it as being the accre- 
tion disc and three years later Osaki (1974) 
suggested that dwarf nova outbursts are trig- 
gered by a disc instability. Few years later 
Hoshi (1979) was the first to find the origin 
of this instability' but it was the talk by Jim 
Pringle at the Sixth North American Workshop 
on Cataclysmic Variables, held at Santa Cruz 
in 1981 (unpublished, but see Bath & Pringle 
1982) that stimulated several authors (Meyer & 
Meyer-Hofmeister 1981; Cannizzo et al. 1982; 
Faulkner et al. 1983; Mineshige & Osaki 1983; 
Smak 1984a) to publish articles with elaborate 
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' His paper singled out the convection as the main 
agent of the instability, while today one would rather 
attribute it to the change in opacities, the two being 
obviously related 



models of unstable dwarf nova discs. Soon af- 
ter Smak (1984b) published a review article on 
dwarf-nova outbursts which still contains al- 
most everything that one needs to know about 
the subject. The more recent review article by 
Lasota (2001) can be considered as just an ex- 
tended commentary to Smak's 1984 paper; the 
only major difference being that it addresses 
also the model of outbursting discs in X-ray bi- 
naries containing black holes or neutron stars. 

At present this 30 year old model, the 
Disc Instability Model (hereafter DIM) is gen- 
erally believed to provide a correct descrip- 
tion of dwarf novae and X-ray transients. The 
word "believe" is used here in its basic sense, 
i.e. "accept that (something) is true, especially 
without proof (OED). Indeed, many counter- 
example to this believe are ignored by most of 
the interested community. Or one should rather 
say: by the community that should be inter- 
ested in these matters. 

The DIM in its standard, or rather original, 
version assumes that, during the outburst cycle, 
the mass transfer from the secondary star, i.e. 
the rate at which the disc is fed with matter, is 
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constant. This simplifying assumption had also 
the advantage that it was clearly demarking the 
DIM from its original rival: the Mass Transfer 
Instability (MTI) model that has been finally 
discarded, the mechanism able to trigger it 
in the stellar companion remaining elusive. 
However, the assumption about the constancy 
of mass-transfer rate encounters two basic dif- 
ficulties. First, the model with this assump- 
tion produces regular, periodic outburst pat- 
terns that have never been observed. Second, 
real CVs show fluctuations, sometimes very 
large, of the rate at which matter is transferred 
from the companion star. It is clear therefore 
that the constant M assumption must be aban- 
doned but this adds a free function of time to 
the model. 

In addition the DIM is not able to repro- 
duce the so-called superoutbursts, observed in 
SU UMa stars, whose amphtude is larger and 
duration longer than observed in normal out- 
bursts. Some systems, such as the celebrated 
short-period binary WZ Sge show superout- 
bursts only. To account for superoutbursts the 
DIM must be substantially modified. Two main 
possibilities have been explored. The first con- 
sists in adding an additional source of vis- 
cosity: tidal torques acting on the supposedly 
deformed eccentric disc. An important evi- 
dence in favour of this eccentricity is sup- 
posed to be the presence of so-called super- 
humps modulated with a period slightly dif- 
ferent from the orbital. The resulting increased 
viscosity causes a prolonged accretion episode 
resulting in large amplitudes and long duration. 
This Tidal Thermal Instability model (TTI; 
see Osaki 2005, and references therein) has 
encountered several serious difficulties when 
compared with observations (see e.g. Buat- 
Menard & Hameury 2002; Hameury & Lasota 
2005, 2006) but probably the strongest chal- 
lenge to its application to superoutbursts has 
been the paper by Smak (2009) in which he 
showed that the presumed evidence of disc ec- 
centricity during superoutbursts resulted "ei- 
ther from errors, or from arbitrary, incorrect as- 
sumptions". The second possibiUty of includ- 
ing superoutbursts into the DIM framework 
consists in explaining their amplitudes by an 
enhancement of the mass-transfer rate during 
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Fig. 1. Response of a steady accretion disc to a 
mass transfer rate variation by a factor 20. The pa- 
rameters used are those of EG Cnc, Irradiation is 
included as is the presence of a hole in the disc as 
well as the existence of a tidal instability. Initially 
M2 = 1-5 X lO'^gs"; at t - 10; M2 increases to 
3.0 X lO'^gs'', and at t = 30 returns to its quies- 
cent value. The top panel shows the mass accretion 
rate onto the white dwarf (solid line) and the mass 
transfer rate from the secondary (dashed line); the 
other panels show the outer disc radius, its mass and 
visual magnitude. (Buat-Menard & Hameury 2002) 



the initial phases of an outburst that otherwise 
would have belonged to the normal category. 
This model has been also proposed by Osaki 
(1985) who, however, has quickly abandoned 
it in favour of the TTI model. Hameury et 
al. (1997) apphed the enhanced mass-transfer 
(EMT) model to WZ Sge whose fluence and 
very long (30 years) recurrence time (after the 
detection of its first outburst it was considered 
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to be a nova) could be reconciled with the DIM 
(also in its TTI version) only if the viscosity 
parameter a in quiescence was assumed to be 
one or two orders of magnitude lower than in 
other dwarf novae. Not only the very long re- 
currence could be accounted for only by low- 
ering the value of a but also the mass accreted 
during the outburst could be explained only 
with very low values of the viscosity param- 
eter, the mass-transfer rate being too low to al- 
low accumulating the required amount of mass 
despite the exceptionally long recurrence (ac- 
cumulation) time (Smak 1993). Buat-Menard 
& Hameury (2002) combined the TTI and EMI 
models to reproduce the lightcurve of EG Cnc 
(Fig. 1). 

Smak (2008) proposed what he calls "a 
purely observational scenario": Superoutbursts 
are due and begin with a major enhancement 
in the mass transfer rate. During the "Hat-top" 
part of the superoutburst the mass transfer rate 
decreases slowly, causing the observed lumi- 
nosity to decline. The superoutburst ends when 
the mass transfer rate decreases below its criti- 
cal value, resulting in a transition to the quies- 
cent state of the dwarf nova cycle. The cause 
for this major enhancement is supposed to be 
irradiation of the secondary star by the ac- 
cretion generated flux. If the regions close to 
the critical Li point are screened by the disc 
this mechanism is unlikely to work (Viallet 
& Hameury 2008). On the other hand, Smak 
(2011) found in five dwarf novae undergoing 
superoutbursts a direct evidence that irradia- 
tion of the companion is modulated at the su- 
perhump period thus supporting his interpreta- 
tion of these features as being due to irradiation 
controlled mass transfer rate resulting in mod- 
ulated dissipation of the kinetic energy of the 
stream. It seems therefore that irradiation con- 
trols mass-transfer during superoutbursts. One 
way solving the contradiction with Viallet & 
Hameury (2008) is to consider the possibility 
of a non-planar, warped disc. 

2. A short course in the DIM 

In a stationary (constant accretion rate), 
Keplerian accretion disc the effective temper- 
ature falls as R^^'^. This follows from angular 



momentum and mass conservation only (see 
e.g. Frank et al. 2002). The effective temper- 
ature at the iimer disc edge can be written as 

T,,^6xlO'{-^T x-'"K, (1) 
Uo.i^i j 

where L38 is the luminosity in units of 10^^ 
ergs"\ T] = 0.1 770.1 is the accretion efficiency. 
Ml accreting body mass in solar units and 
X = c^R/2GM is the radius measured in units 
of the Schwarzschild radius. Jin ~ 4 x lO'^K 
for a cataclysmic variable (L38 ~ 10"^, 770.1 ~ 
1.5 X lO-\ X ~ 3.3 X 10^) and ~ lO^K for 
a neutron star or a stellar-mass black hole ac- 
creting close to the Eddington limit. In general 
discs around compact objects are hot and fully 
ionized in their inner parts but if large enough 
they will reach temperatures at which hydro- 
gen (or another dominant element) starts re- 
combining. The recombination process causes 
drastic changes in opacities modifying the disc 
cooling mechanism. As a result the disc be- 
comes thermally unstable. This is the essence 
of the DIM: a large enough hot stationary disc 
is always unstable. For a given mass the criti- 
cal size for instabiUty depends on the accretion 
rate only: 

where MEdd = i'Edd/O.lc^ is the Eddington 
accretion rate and Rs = 2GM/c^ the 
Schwarzschild radius. For a given fraction 
of the Eddington accretion rate the instabil- 
ity radius descend deeper into the accreting 
body potential well. This is the reason why, 
even when unstable, accretion discs around > 
10** Mq black holes do not produce dwarf-nova 
type outbursts (Hameury et al. 2009) We will 
illustrate the principle of the DIM by consid- 
ering a sequence of models of discs around 
white dwarfs. In Fig. 2 we represent station- 
ary disc surface-density 2(7?) profiles for three 
values of the central white-dwarf mass 0.4, 0.6 
and 1 .OMq and five values of the accretion rate 
in the units of g/s: log M = 13, 14, 15, 16 and 
17. The viscosity parameter is a = 0.01. Each 
line labeled by a value of the accretion rate 
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Fig. 2. Surface density profiles for equilib- 

rium discs around white dwarfs for three mass val- 
ues: 0.4, 0.6 and 1.0 Mq 

represents a stationary disc solution. The in- 
ner "cut-ofF' is caused by the no-stress inner 
boundary condition vl. - and the number at- 
tached is the values of the white-dwarf radius 
since the disc is assumed to extend down to the 
star's surface. Further away from the edge the 
column density falls like R^^^* (see e.g. Frank 
et al. 2002) until it reaches a critical value 
marked as Y.^m where the slope changes and 
becomes positive. In turn, this positive-slope 
branch encounters a second critical surface- 
density Emax where the slope changes to neg- 
ative. The more familiar, local picture of disc 
equilibria forming the celebrated S-curve can 
be recovered by considering all possible sta- 
tionary solutions at a given radius. Taking a ra- 
dius well inside the discs it is easy to see that 
the M(S) points form an S-curve. In Fig. 2 go- 
ing up from low E and M one encounters first 
points on the thick-line segments of the equi- 
librium curves - they represent a growing seg- 
ment of the M(S) curve (the "lower branch" of 



logRiolcm] 



Fig. 3. Heating and cooling fronts in a helium dom- 
inated Y = 0.98, Z = 0.02 accretion disc. (Figure 
courtesy of Iwona Kotko). 

the S-curve), then one gets to the dotted seg- 
ments whose points will form a decreasing seg- 
ment of the M(2) curve (the "middle branch") 
and finally points collected on the thin contin- 
uous lines will form what is known as the up- 
per branch of the S-curve. Therefore negative 
slopes in E(7?) correspond to stable solutions 
while a solution with surface density increas- 
ing with radius is unstable. At given radius this 
corresponds to the familiar stability/instability 
dM/dl. > 0/dM/dI. < conditions. 

Hence for a given mass-transfer rate the 
stability properties of an accretion disc depend 
on it size: if its outer radius Ro > Rem no sta- 
ble equilibrium configuration exists^. In such 
cases the disc goes into an outburst cycle in a 
hopeless quest for a stable solution. For a white 
dwarf primary this results in a dwarf-nova out- 
burst cycle. 

The disc is brought up to a hot bright 
state by a heating front propagating out- or in- 
wards depending on the mass-transfer rate (in 
Fig. 3 the heating front propagates outwards). 
This phase represents the rise to maximum. 
Propagating through the disc, the heating front 
redistributes^ the temperature (flat profile at 
quiescence) and surface-density (at quiescence 

^ Which has not prevented some astronomers 
from attempting to fit stationary disc spectra assum- 
ing effective temperatures spanning from 30 000 to 
4000 K. 

' The front propagation is a thermal and viscous 
process. 
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increasing with radius) distributions so that at 
peak luminosity the disc has quasi- stationary 
structure with M(R) ~ const, and T^ff ~ 2 ~ 
R"^^'^. Since at maximum M w M^^.J^Ro) - 
the upper critical accretion rate - is by con- 
struction (the disc is unstable) higher than the 
mass-transfer rate, the disc must empty out. 
This occurs on a viscous time-scale: ~ R^ jv ~ 
R^/acsH (using the Shakura & Sunyaev 1973, 
prescription for the kinematic viscosity coeffi- 
cient V = aCsH, with H the discs's semi-scale- 
height). Thus by observing dwarf-nova decay 
from maximum one can (measuring the disc's 
size and estimating its temperature - both pro- 
cedures rather straightforward, especially in 
ecUpsing systems) determine the viscosity pa- 
rameter a. Analyzing dwarf nova light-curves 
collected by van Paradijs (1983), Smak (1999) 
has determined that a » 0.2 in hot, ionized ac- 
cretion discs. It is worth stressing that MHD 
simulations obtain, at best, values an order of 
magnitude lower (see e.g. Sorathia et al. 201 1). 
The ball is clearly in the simulators court. 

The viscous decay from maximum pro- 
ceeds through a sequence of quasi-stationary 
configurations. This pushes the outer disc re- 
gions into the unstable regime. As a result 
a cooling front forms bringing the outer disc 
into the cold quiescent state. During the cool- 
ing front propagation one recovers (for obvious 
reasons) a global disc structure similar to that 
when discussing above the stability of equilib- 
rium discs (compare Fig. 3 with Fig. 2). The 
cooling front brings the disc to a cold state 
(onto the "cold branch" of the S-curve) but 
(contrary to what is sometimes asserted) does 
not govern the decay from maximum. The rea- 
son being that the speed of the cooling front 
is (slightly) smaller than the viscous speed (~ 
vIR) in the hot discs it penetrates into (Menou 
etal. 1999). 

It can be considered that the DIM is rather 
successful in describing the rise to and the de- 
cay from maximum. The same cannot be said 
of the quiescent phase of the dwarf nova out- 
burst cycle (see e.g. Lasota 2001, for a discus- 
sion). 



2.1. SS Cygni is not a dwarf nova (but 
looks like one) 

The cataclysmic variable star SS Cygni is 
the brightest and best observed dwarf-nova. 
Although formally not the prototype of its sub- 
class (it is of U Geminorum type) it is in fact 
the epitome of the dwarf nova undergoing nor- 
mal outbursts. The reproduction by the DIM 
of the two main types of normal outbursts ob- 
served in SS Cyg ("A" and "B" or "long" and 
"short") might require small modulations of 
the mass-transfer rate but the multi-wavelength 
long-term hghtcurve of SS Cyg is quite well 
reproduce by the standard DIM (Schreiber et 
al. 2003) if the disc is truncated and the dis- 
tance is not larger than ~ 100 pc. However, 
the HST/FGS parallax distance to SS Cyg has 
been determined to be 166 ± 12 pc, a result con- 
firmed by the modeling of this system's prop- 
erties by Bitner et al. (2007) whose results are 
consistent with a distance of ~ 140 - 170 pc. 
Schreiber & Lasota (2007) showed that at such 
distance both the accretion rate during outburst 
and the mean mass-transfer rate is too high 
to be compatible with the DIM. With such a 
mass-transfer rate, independent of the DIM, SS 
Cyg should be a nova-like star and show no 
outbursts. In other words: at 166 pc SS Cyg is 
not a dwarf nova. 

Recently, Smak (2010) has challenged part 
of this conclusion by deriving the mass- 
transfer rate from the hot-spot (where the ac- 
cretion stream from the companion encounters 
the outer disc rim) luminosity. His result is 150 
lower than that of Schreiber & Lasota (2007) 
and satisfies the DIM criterion. However, since 
in this new framework the accretion rate in out- 
burst is still the one determined by Schreiber 
& Lasota (2007), explaining the outburst by 
the DIM require enhancing the mass-transfer 
rate during outburst by the same factor of 
150. Smak (2010) concludes that: "Nothing is 
wrong with SS Cyg, nor with theory of dwarf- 
nova outbursts". Although it is difficult to dis- 
agree with the first part of this statement, the 
truth of second is less obvious. Indeed, enhanc- 
ing the mass-transfer rate by a factor 150 is 
a highly non-trivial task. Explaining superout- 
bursts of SU UMa and WZ Sge stars requires a 
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factor 100 "only" (Hameury et al. 1997; Smak 
2004) but outbursts of SS Cyg do not look at 
all like superoutbursts. As mentioned above, 
the mechanism of the putative enhancements 
is rather uncertain but e.g. Smak (2004b) esti- 
mated that irradiation will not work for orbital 
periods longer that ~ 6hr which is rather un- 
fortunate since it is just the orbital period of 
SS Cyg (6.6 hr). 

Smak (2010) explains the visual magnitude 
by a mass-transfer rate M « 6.3 x lO'^'gs"'. 
But this is only a small fraction (0.047+0.005) 
of the total quiescent visual flux. According to 
Bitner et al. (2007) a fraction fd = 0.535 ± 
0.075 of this is emitted by the accretion discs. 
In the framework of the DIM the properties 
of this (non- stationary) disc are rather well 
constrained. Its effective temperature cannot 
be higher than the lower critical temperature 
r^jj = 5210/?y'^ 'M0"4 (Lasota et al. 2008) 
and the (increasing with radius) accretion rate 
should be everywhere lower than 

M„i, = 2.64 X 10'^ a^'^^R]-^^Mf-^\ (3) 

The disc radius R\q is in units of lO''^ cm. 

Let us assume for simplicity that half of the 
disc visual flux of V = 5.6 is emitted by the 
disc in quiescence. This corresponds to a lumi- 
nosity in V: 

Lv« 1.0x10^3 erg s"^ (4) 

One can estimate the disc effective tempera- 
ture corresponding to this luminosity by as- 
suming that it is roughly constant with ra- 
dius as shown by observations of eclipsing 
dwarf-novae and required by the DIM. Using 
the system parameters determined by Bitner 
et al. (2007) one gets for the outer disc ra- 
dius {Rd = 0.9 /?i) 4.7 - 5.5 X 10i°cm (see 
Schreiber & Lasota 2007). R\ is the distance to 
the Li point. Therefore the effective tempera- 
ture is Teff aj 6540 K. This is of course higher 
than the critical temperature. With such a tem- 
perature the disc in SS Cyg cannot be in a quiet 
cold state as required by the DIM. Let us note 
that the value we derived is certainly a lower 
limit, since radii of quiescent discs in dwarf- 
novae are smaller than 0.97? i (see e.g. Harrop- 
Alin & Warner 1996) and we took into account 
only the V luminosity. 



Using more refined methods Smak (pri- 
vate communication) showed that one obtains 
Teff's lower than critical by taking fj+i < 0.45 
(where fd+\ is the luminosity fraction emitted 
by the disc plus the primary). In view of the 
uncertain value of this parameter (Bitner et al. 
2007) this is not an unacceptable modification 
but (because of flux conservation) it implies a 
higher luminosity of the secondary making it 
K2 or earlier in contradiction with the K4.5 in 
Bitner et al. (2007). In other words, at 166 pc 
there is too much flux in SS Cyg to accom- 
modate both the DIM and the observed prop- 
erties of this binary's components. Clearly the 
distance-to-SS Cyg problem requires serious 
investigations, after all it is the prototype dwarf 
nova, but nobody (except for Smak and the 
present author) seem to be interested. So much 
for the Golden Age of Cataclysmic Variables. 

3. Outbursts of the hyperluminous 
X-ray source HLX-1 in ESO 243-49 

Observing a Fast Rising Exponentially 
Decaying (FRED) lightcurve of a newly dis- 
covered X-source one is immediately tempted 
to attribute it to the DIM. This was the case 
of HLX-1 in the halo of the edge-on SOa-type 
galaxy ESO 243-49 which is the brightest 
known ultraluminous X-ray source (ULX; see 
Roberts 2007, for a review) with a maximum 
luminosity >10^^ erg s"' (Farrell et al. 2009; 
Godet et al. 2009). This has not only the 
highest (by an order of magnitude) luminosity 
ULX but is also unique in showing clear, large 
amplitude outbursts (see Fig. 4; and Godet 
et al. 2011) whose spectral behaviour is very 
similar to that observed in steUar-mass Black 
Hole X-ray Binary transients. The luminosity 
of HLX-1 (~ 10"*^ erg s"') and the variability 
timescale (~ 10^ s) imply a huge accretion rate 
onto the compact object, of the order of 10~* 
Mq yr~'. Sustaining such a high mass transfer 
rate excludes wind accretion and implies the 
presence of a stellar companion filling its 
Roche-lobe (if only during part of its orbit) 
losing matter that forms an accretion disc 
around the black hole. With such an assump- 
tion one encounters a fundamental difficulty 
because for a black hole mass assumed to be 
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Fig. 4. The lightcurve of HLX- 1 . The red line rep- 
resents a fit with two exponentials with decay times 
T| = 20 ± 4 days, Xi = 130 ± 3 days; the sec- 
ond outburst is longer than the three other ones with 
Ti = 110 ± 4 days, T2 = 135 ± 3 days. The recur- 
rence time is ~ 376 days. (Figure courtesy of Didier 
Barret). 



around the mass ratio q - M,/Mbh, 

where Mbh and M, are respectively the 
black-hole and stellar-companion masses, will 
typically be ^ ^ 10"'* - 10"^ ^ 1. For such 
a small q, matter circularizes very close to 
the donor star and the 2:1 Lindblad resonance 
appears at a radius ^ 0.63a (where a is the 
orbital separation) within the primary's Roche 
lobe (Lin & Papaloizou 1979), affecting the 
accretion disc formation and structure. The 
only thing that is certain is that the standard 
formulae for the sizes of the elements of the 
binary cannot be used in such a case. For lack 
of better solution one assumes the disc forms 
and has a size Ru ~ a. 

As explained in Sect. 2, to be unstable the 
disc must have a radius larger than the critical 
value of Eq. (2), hence for HLX-1 with lO"* Mq 
and M - lOMEdd (Lmux ~ i'Edd) the condition 
is: 



Rr, > R„n = 2.7 X 10'-^ M; 



^/^m;/^cm 



(5) 



where M4 is the mass in units of M© and mio - 
(lOM/MEdd). (The value given by Eq. (5) is 
slightly different from that of Eq. (5) in Lasota 
et al. (201 1) because of different forms of cri- 



teria used - the formulae for critical quantities 
are obtained from fits to equilibrium S -curves 
and often not very exact.) The actual disc can 
be much larger than this limit because: i.) the 
heating front bringing the disc to maximum lu- 
minosity could have stopped before reaching 
the outer edge; ii) in Eq. (5) we have assumed 
a non-irradiated disc whereas the disc in HLX- 
1 is certainly irradiated and therefore stabilized 
(see e.g. Dubus et al. 1999). Hence a larger crit- 
ical radius. 

At first sight there seem to be nothing 
wrong about the radius required by Eq. (5). As 
shown in Lasota et al. (201 1) the Roche geom- 
etry implies a mean density of the companion 
star p X 0.057 Pf g cm^l From Eq. (5) it fol- 
lows that the orbital period is ^ 23 days. The 
secondary star mean density of this hypothet- 
ical companion is lO"'* gcm"^ implying a red 
giant or a massive supergiant - a very reason- 
able option. 

However, even a quick look at the 
timescales completely shatters this optimistic 
conclusion. Indeed, the typical outburst 
timescale will be linked to the viscous 
timescale of the minimum critical radius 



R'- 



t^i, = :_ a; 115 a-^T^^R\'^M\'^yeeiTS 



(6) 



T4 is the disc temperature in units of 10'* K (we 
assume for simplicity that it is the temperature 
at Rent)- Therefore the variability timescale is 
much too long for any reasonable set of pa- 
rameters. The variability in HLX-1 cannot be 
related to the processes described by the DIM 
since that it is impossible to reconcile these 
timescales with a disc large enough that its 
temperature at the outer radius is lower than 
« 10'* K. Therefore the accretion disc in HLX- 
1 must be hot and thermally stable. 

It seems that the only serious option left 
is that rejected for dwarf novae: the MTI. But 
that's another story (Lasota et al. 201 1). 



4. Conclusions 

The Disc Instability Model in its pure, original 
form describes no real system. Even the proto- 
type dwarf nova star, SS Cyg can be described 
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only if its inner disc is trancated and the mass- 
transfer vary. If its Hubble parallax gives the 
true distance then it cannot be described by 
the DIM at all. It should rather be nova-like 
star, but obviously it is not. The lack of inter- 
est in this fundamental problem is rather sur- 
prising. Although the DIM has been applied, 
more or less successfully, to outbursts of X-ray 
transient binaries, rather surprisingly it can- 
not be appUed to outbursts of the Intermediate 
Mass Black Hole HLX-1 despite them exhibit- 
ing most of the required properties. 
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